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We study the inverse Compton scattering of the CMB photons off high-energy nonthermal elec- 
trons. We extend the formalism obtained by the previous paper to the case where the electrons have 
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tion, we find the same scaling law for the probability distribution function Pi,k{s) as Pi(s) which 
corresponds to the zero bulk motions, where the peak height and peak position depend only on 
the power-index parameter. We solved the rate equation analytically. It is found that the spectral 
intensity function also has the same scaling law. The effect of the bulk motions to the spectral 
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the X-ray and gamma-ray emission models from various astrophysical objects with non-zero bulk 
motions such as radio galaxies and astrophysical jets. 
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I. INTRODUCTION 

The inverse Compton scattering is one of the most 
fundamental reactions which have a variety of applica- 
tions to astrophysics and cosmology. They are, for exam- 
ple, the Sunyaev-Zeldovich (SZ) effects}!!] for clusters of 
galaxies (CG), cosmic-ray emission from radio galaxiesQ 
and clusters of galaxies [3], and radio to gamma-ray emis- 
sion from supernova remnants [H, Therefore, theo- 
retical studies on the inverse Compton scattering have 
been done quite extensively for the last forty years, start- 
ing from the works by Jones 6j, and Blumenthal and 
Gould [7J to the recent works, for example, by Fargion|8j, 
Colafrancescojij [Toj . and PetrukfTTj. 

In particular, remarkable progress has been made in 
theoretical studies for the SZ effects for CG. Wright [12f 
and Rephaeli[l3[ calculated the photon frequency redis- 
tribution function in the electron rest frame, which is 
called as the radiative transfer method. On the other 
hand, Challinor and Lasenbvflij and Itoh, Kohyama, and 
NozawaJIH solved the relativistically covariant Boltz- 
mann collisional equation for the photon distribution 
function, which is called the covariant formalism. Al- 
though the two are very different approaches, the ob- 
tained results for the SZ effect agreed extremely well. 
This has been a longstanding puzzle in the field of the rel- 
ativistic study of the SZ effect for the last ten years. Very 
recently, however, Nozawa and Kohyama 16] showed that 
the two formalisms were indeed mathematically equiva- 
lent in the approximation of the Thomson limit. This 
explained the reason why the two different approaches 
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produced the same results for the SZ effect even in the 
rclativistic energies for electrons. With the formalism, 
studies on various formal solutions and numerical solu- 
tions were presented [l7| . 

Furthermore, Nozawa, Kohyama and Itoh[l8[ (denoted 
paper I hereafter) ap plie d the formalism obtained by 
Nozawa and Kohyama[16j to the inverse Compton scat- 
tering of the CMB photons off high-energy nonthermal 
electrons. This extension is particularly interesting for 
the analysis of X-ray and gamma-ray emissions, for exam- 
ple, from radio galaxies Q and supernova remnants [U, @, 
where the inverse Compton scattering of the CMB pho- 
tons off nonthermal high-energy electrons plays an essen- 
tial role. In paper I, a universal scaling law was shown for 
the redistribution function P\(s) and the spectral inten- 
sity function I{x) under a specific condition for the elec- 
tron distribution which is typically realized. It was shown 
that the spectral intensity function for different energy 
scales (, for example, keV, MeV and GeV) were described 
by one equation with a scaling variable X — a;/47^ lin , 
where x is the photon energy in units of the thermal en- 
ergy of the CMB and 7 m i„ is the minimum value of the 
Lorentz factor of the electron power-law distribution. 

In the present paper, we extend the formalism ob- 
tained by paper I to the electrons of the astrophysical 
objects with bulk motions. As for the CG, the effect of 
the bulk motions was originally obtained by Sunyaev and 
Zeldovichpjj], which is known as the kinematical SZ ef- 
fect. The relativistic corrections to the kinematical SZ 
effect was presented by Nozawa, Itoh and Kohyama[20j. 
We will explore the effect of the bulk motions of the as- 
trophysical objects (, for example, the peculiar velocity 
of the CG) for the high-energy inverse Compton scatter- 
ing. This extension will be particularly interesting for the 
analysis of X-ray and gamma-ray emissions, for example, 
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from radio galaxies with non-zero peculiar velocities and 
various astrophysical jets. 

On the other hand, it is well known that the Solar Sys- 
tem (, i.e., the observer) has a bulk motion with respect 
to the CMB frame. Assuming that the CMB dipole is 
fully motion-induced, we deduce that the Solar System is 
moving with a velocity f3g = vg/c= 1.241 xlO -3 towards 
the direction (£,b) = (264.14° ±0.15°, 48.26° ±0.15°) 0- 
[23j . Chluba et al.[24| calculated the corrections to the 
SZ effect for the CG arising from the bulk motion of 
the Solar System. Nozawa, Itoh and Kohyama[2_5| cal- 
culated the effect in more general way with the Lorentz 
covariant formalism. In the present paper, we apply the 
method developed by Nozawa, Itoh and Kohyama[25] to 
the high-energy inverse Compton scattering and calcu- 
late the effect of the bulk motion of the observer. Thus, 
the effect of the bulk motions for both the astrophysical 
object and observer will be derived in the present paper. 

Before closing the present section, it should be empha- 
sized the following: In the present approach, we push an- 
alytic techniques as much as possible in order to obtain 
analytic solutions. In contrast to the direct numerical 
calculation, the present approach will have an advantage 
that one may reveal essential physics properties behind 
the numerical results. In the present paper, under a spe- 
cific condition for the electron distribution which is typ- 
ically realized, we will show that a universal scaling law 
is established for the spectral intensity function even if 
one includes the effect of the bulk motions for both the 
astrophysical object and observer. 

The present paper is organized as follows: In Sec. II, 
we derive the analytic expressions for the redistribution 
functions Pk{s, 7) and P\,k{s). Assuming the power-law 
electron distribution, we show that P\,k{s) has the same 
scaling law as Pi(s), where the peak height and peak po- 
sition depend only on the power-index parameter. We 
calculate the rate equation and obtain the analytic ex- 
pression for the spectral intensity function dI(X) / dr. We 
show that dI(X)/dr also has the scaling law, where the 
peak height and peak position depend only on the power- 
index parameter. In Sec. Ill, we derive the the analytic 
expression for the spectral intensity function dI(X)/dr 
which includes the effect of the bulk motion of the ob- 
server. Finally, concluding remarks are given in Sec. IV. 



II. HIGH-ENERGY INVERSE COMPTON 
SCATTERING OFF ELECTRONS WITH BULK 
MOTION 



A. Rate equations in the Thomson approximation 

In Nozawa and Kohyama[16] , it was shown that the co- 
variant formalism [lH and radiative transfer method [l2T| 
were mathematically equivalent in the following (Thom- 



son) approximation: 
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where uj is the photon energy, 7 is the Lorentz factor, and 
f3 and m are the velocity and rest mass of the electron, 
respectively. Throughout this paper, we use the natural 
unit h = c = 1, unless otherwise stated explicitly. For 
the CMB photons, Eq. ([1]) is fully valid from nonrela- 
tivistic electrons to extreme-relativistic electrons of the 
order of TeV region. In paper l(l8|. the high-energy in- 
verse Compton scattering of the CMB photons has been 
studied under the assumption of Eq. JT]). 

In the present paper, we extend the formalism obtained 
in paper I to the electron distribution with a bulk motion. 
Let us suppose that the astrophysical object (, for exam- 
ple, we consider the CG in the present paper) is moving 
with a bulk velocity (3c (=vc I c) with respect to the CMB 
frame. As a reference system, we choose the system that 
is fixed to the CMB in the present section. We discuss 
the effect of the bulk motion of the observer (the Solar 
System) in Sec. III. The z axis is fixed to a line connect- 
ing the observer and the center of mass of the CG. (We 
assume that the observer is fixed to the CMB frame.) In 
the present paper we choose the positive direction of the 
z axis as the direction of the propagation of a photon 
from the observer to the CG. 

The rate equations for the photon distribution func- 
tion n(x) and spectral intensity function I{x) were de- 
rived in Nozawa and Kohyama[16] under the assumption 
of Eq. ([T]). Here, x = u> /ksTcMB is the photon energy 
in units of the thermal energy of the CMB, and s is the 
frequency shift defined by e s — x'/x. We recall the re- 
sults here to make the present paper more self-contained. 
They are given as follows 



sent pap 
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dn(x) 

dr 
dl(x) 

dr 



dsP 1 { S ,(3 c , z )[n(e s x)~n(x)} , (3) 
dsP^s,^) [e-* s I{e s x) - I(x)] , (4) 



dr = n e cjj'dt , 
Pi(s,0 c>z ) = P 1 (s)+p c , z PiM s ) 



(5) 
(6) 



where I(x) — I x 3 n(x), I = (k B TcMB) 3 /27r 2 , n e is the 
electron number density, ut is the Thomson scattering 
cross section, and j3c,z is the bulk velocity of the CG 
parallel to the observer. It should be noted that 0(j3 c ) 
and higher-order contributions were neglected in deriving 
Eq. because (3c,z <C 1 is satisfied for most of the CG. 

In Eq. ([5]), -Pi(s) was calculated in paper I, and P\^k{ s ) 
is the term which appears in the case of non-zero bulk 
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motions. They are denned as follows: 
Pi(s)= f d00 2 1 5 Pe (E)P( s ,0), 

JPrnir, 

PiM»)= f d00 2 j 5 Pe (E)P K (s,0), 

J0 m tn 



P(s,j3) 



r/(Mo,Mo) 



(7) 
(8) 
(9) 



we rewrite Eq. ([10]) as follows: 
P K (s,0)=S(0) \P K (s,0)-P(s,0) 

Pk(s,(3) 



2/3 7 6 



(21) 

3 /(Mo, Mo) . ( 22 ) 



2fr 



/^Mo - /? 2 



Aii(s) (!-/3mo) : 



■/(Mo, Mo) ,(10) 



where P{s,f3) in Eq. J2TJ) is defined by Eq. © and the 
explicit forms were derived in paper I. (Readers may be 
referred to paper I for the explicit forms.) Therefore, the 
main concern in this section is to derive the explicit forms 
for Pk(s,/3). Note that the identity relation 



/(Mo,M ) = o 



1 + MoK, +2( 1 -/W-^o) 



(11) 
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Imp - l 2 = 

(l-/?Mo) 3 7 2 (1 "/3mo) 3 (1-/W 



In Eq. ([T0| . <$(/?) is a factor related to the electron dis- 
tribution function, which is, in general, a function of 0. 
The explicit forms are given by Nozawa and Kohvama(16j 
for three different electron distribution functions. We will 
define the explicit form later in this section. The electron 
distribution function of a momentum p is normalized by 
dpp 2 p e (E)/m 3 = 1. Variables appearing in Eqs. ([7]) 
- Ijlip are summarized as follows: 



(23) 



was used in deriving Eq. (|2ip. It is also important to 
mention that Eq. (|22l) satisfies the following relation: 



P K (s,(3)e- 4s =Pk(s,(3). 



(24) 



In Eq. (f2"2"]l . the integral of Mo can be done analytically. 
One obtains as follows: for s < 0, 



Pk(s,P) 



[A^e* + A 2 {0)e 2s + A 3 (0), 



0-min = (1 - e- |s| )/(l + e-l^) 
M = [l- e s (l-/3Mo)]//?, 

Mi(s) = ' ' 



for s < 



M2OO 



[1 - er s {l + 0)]/0 for s > 
[1 - e~ s (l -0)}/0 for s < 



1 



for s > 



(12) 
(13) 

(14) 
(15) 



32/3 2 7 4 

+ (A/3 + a) (5i(/3)e s + B 2 (0)e 2s + B^e 38 )] , (25) 



and for s > 0, 



Pk(s,I3) = 



[A 3 (0)e s + A 2 (0)e 2s + A x {0)< 



3,s 



The total probabilities for P(s,0) and Pk{s,(3) are 
given by 



32/3 2 7 4 

+ (X p - s) (Bi(/3)e s + B 2 {0)e 2s + Bi(,3)e 3s )] , (26) 

where the coefficients are 
1 



dsP(s,j3) = 1, 



+A/3 



dsP K (s,0) = -<5(/3)/3 2 



(16) 
(17) 



Ai(/8) 
A 2 (/3) 



/3 4 7 2 (l+/3) ; 



(9 + 12/3 - 4/3 2 - 80 3 - 3/3 4 ) , (27) 



4M + 2/3 2 ), 



(28) 



/?V(1 - /3) ; 



(-9 + 12/3 + 4/3 2 - 8/3 3 + 3/3 4 ) ,(29) 



where 



In 



1 + 
1-/3 



(18) 



It should be noted that the following useful relations: 



Bx{0) = 
B 2 (0) = 



(3~/3 2 ) 

/34 7 2 
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/3 4 7 2 



(30) 
(31) 



P(s,f3)e- 3s = P(-s,0) 
P 1 (s)e- 3s = P^-s) 

are valid for P(s,0) and Pi(s). 



(19) 
(20) 



B. Pk(s,/3) for extreme-relativistic electrons 

In this section, we derive the analytic expression of the 
frequency redistribution function Pr-(s,/3) for extreme- 
relativistic electrons. In order to proceed the calculation, 



It is clear that Eqs. (|25|) and (|26|) satisfy the relation of 
Eq. 424]). 

Now let us consider the case for electrons of extreme- 
relativistic energies E (= 7771c 2 ) 3> mc 2 . Thus, 7 > 1 
and « 1 are assumed. Under this approximation, one 
can rewrite Eqs. (|2"5"|) and (|26p . and one finally obtains as 
follows: for s < 0, 



Pk{s,j) 



32 7 4 



3e s 



3s„,2 



— -8^+8^7 



(A, 



2e s 



(32) 
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and for s > 0, 
Pk{s,i) 

where 



32 7 4 



(A, 



A 7 = 21n(2 7 ) . 



3e 

_1_ 


33 


2 7 2 


2e 3s " 




7 2 - 


) 



(33) 



(34) 



and the expression Pk(s, 7) was used instead of Pk(s, /?)■ 
It should be noted that only leading-order terms were 
kept in deriving Eq. (1331) . and the relation of Eq. (f2~4")l 
was applied in deriving Eq. (|3"2"j) . Thus, one has the total 
probabilities 



+x -< 4 / 1 \ 

dsPK(s !7 )= 3+O — j , 

' X-y \ T / 

&PKS, 7 ) = -%)+Oh 
where the expression <5( 7 ) was used instead of 5(/3) 



(35) 
(36) 



C. Scaling law of Pi,k(s) for nonthermal electrons 

In order to proceed calculation for practical applica- 
tions, let us specify the electron distribution function. 
High-energy electrons in the supernova remnants and ac- 
tive galactic nuclei, for example, are most likely nonther- 
mal. It is standard to describe the nonthermal distribu- 
tion in terms of the power-law distribution function of 
three parameters: 



N 7 r 
0, 



Imin < 7 < 7 

elsewhere 



(37) 



where 7 is the Lorentz factor and iV 7 is the normalization 
constant. In Eq. (|37[) . a is the power-index parameter, 
j m i n and j m ax are parameters of minimum and maxi- 
mum values for 7 , respectively. It is known by Nozawa 
and Kohyama[16j that the factor 6(j3) is a constant which 
depends only on the power-index in the case of the 7- 
power distribution of Eq. (|37[) . Therefore, we use the 
expression 5 instead of 6(f3) hereafter. 

Then, using the definition of Eq. (f2Tj) , Eq. flSJ) can be 
reexpressed as follows: 



Pi,k(s) = S Pi,k(s) ~ Pi(s] 



(38) 



where P\(s) is defined by Eq. ([7]) and the explicit forms 
were derived in paper I. In Eq. (|3"5j) . the function Pi i k(s) 
is defined as follows: for s < 0, 



PiA*)= / d 7 p e ( 7 )Px(s, 7 ), (39) 

Jmax(7 mi „,e~ s / 2 /2) 



where Pi,k(s, 7) is given by Eq. (|32|) . and for s > 0, 



ic(7min.e s/2 /2) 



d7Pe(7)-Pk-(*,7). (40) 



where Pk(s, 7) is given by Eq. (J33J) . In deriving Eqs. (|39|) 
and ([4*0]) . /3 pa 1 was assumed, and the phase space fac- 
tor 7 2 was absorbed, for simplicity, by the power-index 
a in p e ( 7 ). Since the explicit form for the electron distri- 
bution function has been fixed, the parameter S is now 
determined as follows uM: 



S = a + 2. 



(41) 



In the case of the power-law distribution of Eq. (|37l) , 
Eqs (f3"9"]l and (|4T)j) can be integrated analytically. The 
explicit forms are given as follows: for — 21n2 7maa; < 
s < -21n2 7mm , 



Pi.k(s) = ^iV 7 



3<r + 23 



cr + 5 [2(ct + 5) 
1 / 3ct + 23 



2<r+5 e (o-+5)s/2 



7mtl V 2(a + 5) 
8e 2s 



(7 + 3 

8e 3s 



cr + 1 

for -21n27 mm < s < 0, 



2 o-+3 e (<T+3)s/2 
2 o-+l e (<T+l)s/2 



2s + 41n2 7 
1 



Jmax 



1 



V+ 1 
Jmax 



(42) 



3_ 

32" 
1 

Z7+5 

1 min 
2s 



a + 5 
1 

fmax 



3ct + 23 
2((7 + 5) 

4 m 



7 -+ 5 

' 1)7 +77 



-2s 

imax 



a + 3 V 7 CT+Ci 7mtx 

\ 'mm [max 



(7 + 1 V7 ff+1 7mtx 

\ imin /max 



(43) 



for < s < 21n2 7min , 



Pi ' k(s) = 32^ 



,3s 



(7 + 5 

1 



Imin 1 max 



3cr + 23 
2((7 + 5) 

4 m 



7 



(T+5 



- 2s 

^Jmin 
^max 



^ + :J 'V7^I„ 3 7mtl 
1 1 



(7+1 



imin jmax 



(44) 
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and for 2\n2j min < s < 2ln2"/ max , 



and for < s < 2 In 2"f„ 



Pi,k(s) = -N, 



' f 3(7 + 23 2CT+ 5 e -(<r+5) S /2 



(7 + 5 {2(a + 5) 



1 (3-f23 2s + 

^ 20-+3 e -(cr+3)s/2 _ 1 



(7 + 3 



fmax 



| 8e * f 2 <r+l e -(<7+l)«/2 



(45) 



It should be noted that the normalization constant is 
given by 



JV 7 = (a - 1) 7 ; 



.(7-1 

min 



(46) 



for the case 7 max — > oo. 

Let us now introduce new functions Pk,c{s, R) and 
P' K IC (s, R) in order to express Eqs. (|4"2"1) - (|4"5"]) in unified 
forms, where R = ^min/lmax- Here, C and IC denote 
the Compton scattering and inverse Compton scattering, 
respectively. First, we define Pk,c(s,R) as follows: for 
-21n27 mi „ < s < 0, 



Pk,c{s,R) 



3fr-l) 
1 - R"- 1 



a 3s 



(7 + 5 IV (7 + 5 



3cr + 23 



4.s 



(1 - P CT+5 ) + 8P CT+5 Inpj 

/| c 2s s 

H — - (l — -R 

(T+ 1 V 



1-P CT+3 

(T + 3 v 

and for < s < 2hi(j max /j min ), 



, (47) 



Pk,c{s,R) 



3(a-l) [ e 3s (3o + 23 _ (ct+5)s / 2 



1 - P° 



, : 5 ( 3a + 23 
cr + 5 



-P 

4e 2s 



(7 + 5 ]_ ct + 5 
- 4s - 81nP 



+ 



_ ( e -(^+3)s/2 _ d(7+3 

(7 + 3 V 

(ct+1)s/2 _ ^ff+l 



(7 + 1 



(48) 



Similarly, P' KIC {s,R) is for -21n(7 maa: /7 min ) < s < 0, 



3(7 + 23 



(7 + 5 (7 + 5 



,(<7+5)s/2 



_ flCT+5 /3a + 23 4s _ 81ni? 

V (7 + 5 



4e 2s / 

(7 + 1 V 



(<r+3)s/2 _ ^(7+3 



t/2 _ ^(7+1 



(49) 



P'kjc{s,R) 



3(t~1) 
1 - 



(7 



/3<7 + 23 



5 tV o- 



4s 



x (1 - i? cr+5 ) 
4e 2s 



1 -P 



8R a+t> \nR 

<7+3\ | 



(7 



(7+1 



1-P° 



.(50) 



It is straightforward to show that 

P K ,c{s 1 R)e- 4s = P' KJC (-s,R) (51) 

is satisfied by Eqs. (l4"71)- ([50")) . 

Comparing Eqs. (|4"2" ]) -(|4"5 j) with Eqs. (|47 ]) -(|5n ]l , the 

probability distribution function P\.k{s) is described as 
follows: 

Pk.ic (s + 2 In 2 7mi „ , R) for s < 

PiA*) = < . (52) 

P^c (a - 2 In 2 7mm , P) for s > 



where 



P K .ic{s,R) 



256 7 .; 
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ravn 



P' KJC (s,R). 



(53) 



Combining the result of Eq. (|52|) for A, if (s) with the 
previous result for Pi(s) in paper I, one finally obtains 
the total contribution P\,k (s) defined by Eq. ([55]). 

Let us now consider the case P = Jmin/lmax *C L We 
fix 7mocc = 10 8 throughout the paper. In Fig. 1(a), we 
plot Pi, x (s) defined by Eq. ([38]) as a function of s for a 
typical value a = 2.5. The solid curve, dash-dotted curve, 
dashed curve, and dotted curve correspond to jmin = 10, 
10 2 , 10 3 , and 10 4 , respectively. It can be seen that the 
height of Pi, if (s) is independent of J m in- I n Fig. 1(b), 
we plot the same curves as a function of new variable sc 
which is defined by 



sc 



21n27„ 



(54) 



In Fig. 1(b) the four curves are totally indistinguishable, 
which exhibits a scaling law for Pi, if (s). The reason for 
this scaling law is as below. For large j m i n 3> 1, as 
shown by Figs. 1(a), 1(b), and Eqs. ([52]) and ([53]). the 
probability distribution function Pi,^-(s) is dominated by 
Pk,c(sc,0): i.e. by the Compton scattering process. 

Before closing this subsection, we study the a de- 
pendences on the peak position s pea k and peak height 
Pi, if (speak)- As shown in Figs. 1(a) and 1(b), the 
dependence of Pi,if (s) is described by Eq. (|5"4"|) . namely, s 
— sc + 2 In 2^ m i n . Therefore, wc define the peak position 
by 



^peak 



sk{<?) + 21n27„ 



(55) 



where Sif(c) depends only on a. The peak position is 
calculated by solving the equation 



dPl,K{») 



= 0. 



(56) 
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FIG. 1: Plotting of Pi,k(s) and Pi,k(s C ) for a = 2.5. Fig- 
ures 1(a) and 1(b) are Pi,k{s) and Pi,k(so), respectively. 
The solid curve, dash-dotted curve, dashed curve, and dotted 
curve correspond to 7 m i„ = 10, 10 2 , 10 3 , and 10 4 , respectively. 



The analytic expressions for sk (c) in the first-order and 
third-order approximations are given as follows: 

(a- l)(cr 2 + 4a + ll) 



SK,Zrd{&) 



((7 - 3)(5(7 2 + 24(7 + 43) ' 

1 



(57) 



25cr + 101 



13cr + 41 



x/(g + 3)A 3 + B 2 — B 

((7 + 3) 2 



y/(a + 3)A 3 + B 2 +B 
(a + 3) 2 



1/3 



1/3- 



(58) 



A = 81(7 3 + 637cr 2 + 2119(7 + 3643 , (59) 
S = 803cr 5 + 18351cr 4 + 148554cr 3 + 590290ct 2 

+1282323(7 + 1318911 . (60) 

We also solved Eq. (f5"6"]) numerically and obtained the nu- 
merical solution SK,num(&)- In Figs. 2(a) and 2(b), we 



FIG. 2: Plotting of sk(ct) and Pi,K(s pe ak)- Figures 2(a) and 
2(b) are s_r-(<7) and Pi,K{s pea k), respectively. The dashed 
curve, dash-dotted curve, and solid curve correspond to the 
first-order approximation, third-order approximation, and nu- 
merical solution, respectively. 



plot srt((t) and Pi,K(speak), respectively. The dashed 
curve, dash-dotted curve, and solid curve correspond 
to s K> i st (a), s K: 3rd(cr) and s K>num (a), respectively. In 
Fig. 2(a), the solid curve, and dash-dotted curve are 
indistinguishable, and three curves are almost indistin- 
guishable in Fig. 2(b). It can be seen from Figs. 2(a) and 
2(b) that the third-order approximation is sufficiently ac- 
curate for the present purposes. 



D. Scaling law for spectral intensity function 

Let us now solve the rate equations of Eqs. ((3]) and (j4| 
with the results for Pi(s) and Pi,k(s)- First, one can 
rewrite Eq. (QJ as follows: 



dl{x) dli{x) , _ dI K {x) 
~ - PC,z- 



Ot 



Ot 

(1 - Pc,zS) 



dr 
dlijx) 
8t 



d~I K {x) 
dr 



(61) 
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where dli{x)/dr is the result obtained in paper I, and 
dlxix) / dr corresponds to non-zero bulk motions. In de- 
riving Eq. (foTj) . we used the following relation: 



dr 



= 5 



di K (x) dhix) 



(h 



(62) 



Therefore, our main concern in this section is to calculate 
dlK(x)/dr, which is defined as follows: 



dI K (x) 
dr 



dsAzW^/ote^-IoW] , (63) 



where I (x) = I x 3 /{e x - 1) and J = (fc B T CM s) 3 /2tt 2 , 
because we consider the CMB photons for the initial dis- 
tribution. It should be noted that the useful relation 
Pi,K(s)e~ 4s = Pi : k{—s) was used in deriving Eq. (|63|) . 
Similarly, Eq. can be rewritten as follows: 



dn(x) 
dr 



n a x\ dn i( x ) , a x dfl K{x) 

(1 - PC,zO) — h Pc.zO ■ 



Ot 



(64) 



where dni(x)/dT is the result obtained in paper I, and 
dn K {x)jdT is defined as follows: 



dn K (x) 
dr 



dsPi, K (s) [n {e s x) - n (x)} , (65) 
I dI K (x) 



IqX 3 dr 



(66) 



where n (x) — l/(e x — 1). 

For the inverse Compton scattering of the CMB pho- 
tons off high-energy electrons, we are interested in high- 
energy spectra such as X-rays (~ keV) and gamma-rays 
(~ MeV). Therefore, one can safely assume 



ksTcMB 



> 1 



(67) 



for scattered photons. For the 7 parameters, we assume 
the same condition used in the scaling law for P\,k (s), 
namely, 



1 < 7m m < 7 II 



(68) 



Under these assumptions, Eqs. (|63l) and (|65j) are much 
simplified, and can be solved analytically. In Ap- 
pendix B of paper I, we have shown the derivation for 
dli(x)/dr and dni(x)/dT in detail. One can also calcu- 
late 8Ik{x) / dr and dfiK{x) / dr in a similar manner. 
The final results are as follows: 



d~I K {X) 
dr 
I 



f 



3(<t - 1)7 



4(a 2 



x 

3ct + 23 
a + 5 

-4ct+11) 



X- 1 - 



dt- 



1 



4 In- 



X 



x e* - 1 
4 t 



1 t 2 



( ( j + l)(a + 3)( ( 7 + 5) 2 X(-- 1 )/ 2 Jo 



a + 3X <t + 1X< 

X £(<r+5)/2 

dt— 



(69) 



dn K (X) _ 1 I dI K (X) 
dr _ 64 7 L„/oA3 dr ' 

where X is the scaling variable defined by 

X = 



4 7 



2 

rain 



(70) 



(71) 



In order to compare the present results for dIi((X) / dr 
and dfiK{X) / dr with dIj(X) / dr and dni(X)/dr, we re- 
call the results of paper I. They are given as follows: 



dli(X) 
dr 



3(a - 1)I 



x 



dt 1 
Te'-l 



a - 1 



(7 + 5 er + 3 \er + 3 
2(cr 2 +4(7 + 11) 



2 In 



t \ t 



1 



((7 + l)((7 + 3) 2 ((7 + 5)X( ff - 1 )/ 2 7 e 



x 7 x ' cr + 1 x 2 

, (721 



X x(a+3)/2 

dt' 

1 



1 1 

64 7 £ in 7 X3 dr 



(73) 



d ni (X) 
dr 



Comparing Eqs. (|69J) and J70J) with Eqs. (JTSj) and (j73|) . 

respectively, one finds as follows: 



<jfo(X) _ 1 dlijX) 

dr x dr ' 

dn K (X) 1 dn/(X) 



(/T 



dr 



(74) 
(75) 



Thus, it is found that dIx{X) / dr and dhx/dr are sup- 
pressed by a factor 1/x compared with dIj(X)/dr and 
dni/dr, respectively, and they are negligible for x 1. 
Finally, one obtains the total contribution as follows: 



dl(X) 

dr 
dn{X) 

dr 



= (1 - Pc.zS) 
= (1 - Pc, z 5) 



dli(X) 
dr ' 

dni{X) 
dr 



(76) 
(77) 



As seen from Eqs. ([75]) and ([77]) . the same scaling- law 
shown in paper I is valid to the present case which in- 
cludes the effect of the bulk motions. Moreover, it should 
be noted that the effect of the bulk motions is approxi- 
mately expressed by the factor —f3c, z S, where /3c, z is the 
bulk velocity parallel to the observer, S — a + 2, and a 
is the power-index of the power-law electron distribution 
function. It turns out that the effect of the bulk motions 
is small for most of the astrophysical objects. For exam- 
ple, fic,z ~ 1/300 for a typical CG and 5 — 6 for a typical 
value of (7 = 4, one has —(3c,z5 ~ —2%. Therefore, the 
astrophysical applications of the scaling laws suggested 
in paper I are still fully valid even for the astrophysical 
objects which have the bulk motions. 
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III. EFFECT OF OBSERVER'S MOTION 
A. Kinematics 

In the present section, let us discuss the effect of the 
observer's motion. It is well known that the Solar Sys- 
tem has a bulk motion with respect to the CMB frame. 
Assuming that the CMB dipole is fully motion-induced, 
we deduce that the Solar System is moving with a ve- 
locity 0s = vs/c — 1.241 xl0~ 3 towards the direction 
(I, b) = (264.14° ± 0.15°, 48.26° ± 0.1 5°) [211-13. Chluba 
et al.[24[ calculated corrections to the SZ effect for the 
CG arising from the bulk motion of the Solar System. 
Nozawa, Itoh and Kohyama[25| calculated the effect in 
more general way with the Lorentz covariant formalism. 
In the present paper, we apply the method developed by 
Nozawa, Itoh and Kohyama 25] to the high-energy in- 
verse Compton scattering and calculate the effect of the 
bulk motion of the observer. (Readers may be referred to 
the paper for the details.) 

Let us suppose that the observer's system (the Solar 
System) is moving with a velocity 0s{= vs/c) with re- 
spect to the CMB. The z axis is fixed to a line connecting 
the observer and the center of mass of the CG as defined 
in Sec. II. Note that variables in the Solar System will be 
denoted by the subscript S, unless otherwise stated ex- 
plicitly. First, the photon energies lo and los are related 
by the Lorentz transformation: 



7sws(l 
1 



0sVs) , 



is 



v/1 



■ft 



(78) 
(79) 



where [is = —0s'ks, and /3s is a unit vector in the 
direction of 0s- The photon unit wave vectors k and ks 
are related by[26( 



throughout the present paper. In this approximation, 
Eqs. ([75]). ([501 and ([Hi]) are much simplified as follows: 



LU = LU S {1 ~ 0SHS) : 

k = k s , 

0c = C + 0s ■ 



(83) 
(84) 
(85) 



Inserting Eqs. ([B"3"]) and ([541 into Eq. ([8"2"]) . one obtains as 
follows: 



n(x) = ns(xs) 



1 



e x s {l-fistis) — I 



1 



e xs - 1 



^^rif+^i), (86) 



where x s = uj s /k B T C MB, and 0sHs = ~Ps,z- In 
Eq. (|86[) . the second term is the first-order correction 
due to the bulk motion of the observer. 



B. Scaling law for spectral intensity function 

Now, our task in the present section is to derive the 
analytic expressions for the spectral intensity function 
and photon number distribution function with variables 
in the observer's system. This can be done by inserting 
Eq. ([55]) into Eqs. ([75]) and {7_7J). One finally obtains 
the following expressions which include the effect of the 
bulk motions for both the astrophysical object (CG) and 
observer (Solar System): 



dl(Xs) dljjXs) dIs(X s 
— j- — = (1 - Pc.7.0) — Ps. r 

dn(X s ) 



dr 



dr 

n r x\ dn i( x s) a 

= (1 - PC,zO) j- PS,; 



dr ' 
dns(Xs) 



(87) 



dr 



k S -0S+ Ps 
1 + k s ■ 0s , 



0s 



ks - (}s ■ 0s)0s 
7s(l + k s ■ 0s) 



(80) 



Similarly, the velocity of the CG in the CMB frame 0c 
is related by the velocity of the CG in the Solar System 
0'c by the following relationship [2 6]: 



0c 



P'c-Ps + Ps]^ , 0'c - {0' c ^s)0s _ m 



1 + 0'c-Ps 



ls {l + 0' c -0s) 



n(u,h) =n s (o;s,fcs). 



(82) 



where X s — xs/47^ in . Note that 0c, z is the bulk veloc- 
ity of the CG in the CMB frame. The explicit forms for 
dIs(Xs)/dr and dns(Xs)/dr are given as follows: 



dlsjXs 
dr 



Thus, the kinematics for the CMB, CG and Solar System 
are well defined. Then, the initial (Planckian) photon 
distribution functions in the CMB frame and the Solar 
System are related as follows t 25j: 



For the practical calculation we now assume a con- 
dition 0s <C 1, which is well satisfied. Therefore we 
keep only the first-order terms in 0s in the calculations 



3(a - 1)7 
1 



dt- 



2(a 2 4 



a - 
4cr - 



3 

11) 



( ( t + 1)( ( j + 3) 2 ( ( 7 + 5)^ ,t - 1)/2 



x s ~>'-l) 2 
t 

1 

7V2 / 



o- - 1 „ n t 

2 In — 

o- + 3 X s 



1 t 2 
+ 1X1 



Xs t (a+5)/2 e t 

dt-rz 

(e* - l) 2 



dns(Xs) 
dr 



1 1 dI s (X s ) 



64 7 ; 



T Y'S 
min 1 0^S 



dr 



,(89) 



(90) 



As seen from Eqs. ([87]) and ([89]) . the spectral intensity 
function has the scaling law as the previous case even 
after including the effect of the bulk motions for both the 
astrophysical object (CG) and observer (Solar System). 
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FIG. 3: Plotting of dI I {X s )/dT and l3 s , z dI s (X s ) /dr for 
Ps,z = 1.241 x 10" 3 as a function of X s . The solid 
curve is dIi{Xs) /dr and the dashed curve corresponds to 
Ps,zdIs(Xs)/dr. In the calculation, a typical value a — 3.5 
was used for an illustrative purpose. 

Therefore, it should be emphasized that the scaling law 
is universal. 

Furthermore, the effect of the bulk motion of the ob- 
server is small compared with the leading-order term 
because of the smallness of (3s- In Fig. 3, we plot 
dI I {X s )/dr and p StZ dI s {X s ) / 'dr for (3 S . Z = 1.241 x lfr 3 
as a function of Xs- In the calculation, a typical value 
(7 = 3.5 was used for an illustrative purpose. It is clear 
from Fig. 3 that the effect is quite small. The correction 
is roughly 0.4% at the peak position. 



IV. CONCLUDING REMARKS 

In paper l[l8j], we studied the high-energy inverse 
Compton scattering of the CMB photons off nonther- 
mal electrons with the formalism derived in the Thomson 
approximation [Iff. As for the nonthermal electron distri- 
bution function, a standard power-law distribution func- 
tion of three parameters were adopted: the power-index 
a, minimum value "f m im and maximum value j m ax of the 
distribution range. For the case 7 m j T[ 3> 1, a scaling law 
in the probability distribution function Pi(s) were found, 



where the peak position depends on s — 2ln2j m i n , and 
the peak height depends only on the power-index param- 
eter. The spectral intensity function Ii(x) was also cal- 
culated. For the case of high-energy photons of x ^> 1, a 
scaling law in dli{x)/dr were found, where the function 
depends on a new variable X — xj (47„ in ). The P ea k P°~ 
sition and peak height depend only on the power-index 
parameter. The 7 m i„ dependence of dli(X) / dr is in- 
cluded in the variable X . The formalism was applied to 
the observations of the spectral intensity function in the 
X-ray and gamma-ray energy regions. It was found that 
the observations in the X-ray and gamma-ray regions had 
sensitivities of 7mm=500 ~ 3xl0 3 and 7 m i„=16xl0 3 ~ 
95xl0 3 , respectively. 

In the present paper, we have extended the formalism 
obtained in paper I to the cases where the astrophys- 
ical objects have the bulk motions with respect to the 
CMB frame (, for example, the peculiar velocity of the 
CG). The extension will be particularly interesting for 
the analysis of X-ray and gamma-ray emissions, for ex- 
ample, from radio galaxies with non-zero peculiar veloc- 
ities and various astrophysical jets. 

First, we have derived the analytic expressions for the 
redistribution functions Pfr(s, 7) and Pi.k(s). Assuming 
the power-law electron distribution, we have shown that 
Pi.k(s) also has the same scaling law as Pi(s), where 
the peak height and peak position depend only on the 
power-index parameter. 

Then, we have calculated the rate equations and 
obtained the analytic expressions for dI(X)/dr and 
dn(X)/dr which include the effect of the bulk motions. 
It has been found that the same scaling-law shown in 
paper I is valid to the present case which includes the ef- 
fect of the bulk motions, where the peak height and peak 
position depend only on the power-index parameter. 

It has been found that the effect of the bulk motions 
is approximately expressed by the factor —(3c,z3, where 
f3c,z is the bulk velocity parallel to the observer and 
5 = a + 2 is a factor related to the electron distribution 
function. It has been shown that the effect of the bulk 
motions is small for most of the astrophysical objects. 
For example, (3c, z ~ 1/300 for a typical CG and (5 = 6 
for a typical value of <r = 4, one has —f3c,z& ~ — 2%. 

We have also calculated the effect of the bulk motion 
of the Solar System by applying the method developed 
by Nozawa, Itoh, and Kohyama[25J. The analytic expres- 
sion for the spectral intensity function and photon num- 
ber distribution function have been derived. It has been 
found that the same scaling laws is valid for the spec- 
tral intensity function. Furthermore, it has been found 
that the effect of the observer's motion is small compared 
with the leading-order term because of the smallness of 
the velocity (3s- The correction is roughly 0.4% at the 
peak position for a typical value of a = 3.5. 

In conclusion, the astrophysical applications of the 
scaling laws suggested in paper I are still fully valid even 
if the astrophysical object and observer have the bulk 
motions. 
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